MAGNETIC MONOPOLES, ELECTRIC NEUTRALITY 
AND THE STATIC MAXWELL-DIRAC EQUATIONS 



CHRIS RADFORD AND HILARY BOOTH 

Abstract. We study the full Maxwell-Dirac equations: Dirac field with 
minimally coupled electromagnetic field and Maxwell field with Dirac current 
as source. Our particular interest is the static case in which the Dirac current 
is purely time-like - the "electron" is at rest in some Lorentz frame. In this 
case we prove two theorems under rather general assumptions. Firstly, that if 
the system is also stationary (time independent in some gauge) then the 
system as a whole must have vanishing total charge, i.e. it must be 
electrically neutral. In fact, the theorem only requires that the system be 
asymptotically stationary and static. Secondly, we show, in the axially 
symmetric case, that if there are external Coulomb fields then these must 
necessarily be magnetically charged — all Coulomb external sources are 
electrically charged magnetic monopolcs. 
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I. Introduction 

The Maxwell-Dirac equations are the classical field (or more "traditionally" , first 
quantised) equations for electronic matter. Historically, only the linearised 
equations (where the Dirac current is ignored as a source for the Maxwell 
equations) have been studied in detail - for a comprehensive survey of the Dirac 
equation with various potentials see Thaller [|l| . The lack of past interest in the 
full Maxwell-Dirac equations is partly due to the very difficult nonlincarities of 
the equations. More importantly , the classical problem was swamped by the 
extraordinary success of QED. 

The difficult nature of these nonlinear equations has meant that the existence 
theory has only recently been enunciated - some highlights in this development 
might be Gross Q, Chadam 13], Georgiev Q, Esteban et al ||, and Bournaveas 
Q. This work culminated in a tour de force of nonlinear functional analysis, the 
global existence proof of Flato, Simon and Taflin [Q. 

Our aim in studying the Maxwell-Dirac system is to look for possible non-linear 
behaviour which would not be apparent in perturbation expansions. The 
particular solutions found in || and [[Io| exhibit just this sort of behaviour - 
localisation and charge screening. See also Das || and the recent work of Finster, 
Smoller and Yau jn|. 

The static Maxwell-Dirac equations were first written down in ||. In the present 
work we use this formulation to prove two theorems. Firstly, that the stationary, 
static Maxwell-Dirac system must have vanishing total charge, this is done in 
section IV. The second theorem proves that, in the axially symmetric case an 
external Coulomb field must have an associated magnetic charge - external 
Coulomb fields must be electrically charged magnetic monopoles. This theorem is 
proved in section VI. 

II. The Static Maxwell-Dirac Equations 

In standard notation the Dirac-Maxwell equations are 

J a (da — ieA a )ip + imip = 

Fa/3 = Afj, a — A a fi 

d a F a f3 = —4irejp = —Aireipjpip 

In H the 2-spinor form of the Dirac equations was employed to solve for the 
electromagnetic potential, under the non-degeneracy condition j a j a =/= 0. 
Requiring A a to be a real four- vector then gave a set of partial differential 
equations in the Dirac field alone, the reality conditions. 
For 2-spinors u A and v B we have 

i> = y !J| ) , with 

ucv c 7^ (non-degeneracy). 
The electromagnetic potential, 

A AA = -l—L A d BA u B + u A d BA v B + ^(u A u A + v A v A ) 



e(u c v c ) \ " ^2 

The reality conditions, 

d AA (u A u A ) = - ^(^c-^c) 
d AA {v A v A ) = ^L(u c v c ~u c v d ) 
u A d AA v A ~v A d AA u A = 
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The Maxwell equations, 

d a F a (3 = -Awejp = -Aire<Tp A {u A u A + v A v A ). 

These equations constitute the Maxwell-Dirac system. 
We next impose the static condition. 

Definition 1. The Maxwell-Dirac system is said to be static if there exists a local 
Lorentz frame in which the Dirac current vector is purely time-like, i.e. j a — j°5q . 

Imposing this condition one quickly finds that 



ix V2a 



OAA- 



' A' 



with x a re al function. 



The current vector is now 

the gauge is fixed by the choice, 



As noted in 



(1) 

with X, Y, and n real functions on IR 
Defining the null vector L, 

1 

71 



u° = Xe^ x + rt) 

4 



L = (a a . A u A u A ) = (L 



V), with V 



V2 



(X 2 



Y J ) and 



V 



(2XYcosn, 2XYsinn,X 2 -Y 2 



our equations become, 



!-(* 2 
dV 



v.v 

dV 
~dt 

With electromagnetic potential 



Y A ) = 

2m(X 2 + y 2 )sinx 

+ (v x )xy = o. 



^0 m 

A = — cos x 

e 

1 



(X 2 -Y 2 ) On 



(v x ).v 



2e{X 2 + Y 2 ) 



2e(X 2 

'dx 
dt 



+ Y 2 ) dt 2e(X 2 + Y 2 ) 
V + (X 2 - Y 2 )Vn - V X V 



where A = (A 1 , A 2 , A 3 ). 

The full system is given by the above two sets of equations and the Maxwell 
equations. 

One further condition we want to impose is that of stationarity. 

Definition 2. The Maxwell-Dirac system will be called stationary if there exists a 
gauge in which ip = e lu)t <j), with the bi-spinor <j) independent oft. Such a gauge will 
be referred to as a stationary gauge. 

We now have the following simple lemma. 

Lemma 1. The static Maxwell-Dirac system is stationary if and only if in the 
gauge given in (|]) ; 

— = and — - = 0. 

at at 



In the stationary case we also have, 

dx 
dt 



and 



dV 

dt 
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Proof. If the system is stationary there exists a gauge transformation such that 
Consequently, 



u A -> e li u A = e luJt C A , with ( A independent of t. 



y e i(x->7) _ e i(a>t-£) /-l _ 
So, |X| = |C°| and |F| = If 1 ], both independent of t. We also have, 

at"dt = 2( "-^ } - 

So that, ^ = 0. The argument is easily reversed to get the converse 

statement. □ 



III. Isolated Systems 

An isolated system is one for which all sources are contained in some ball B p 
(p < oo) and for which the fields die off as \x\ — r — > oo. 

In what follows we will be considering stationary Maxwell-Dirac systems. For 
such systems we would expect, in an appropriate stationary gauge, that A a 
should be 0(1/ r) as r — > oo. We will also need to impose some decay conditions 
on ip as r — > oo in order to appropriately define an isolated Maxwell-Dirac system. 
The best language for the discussion of such decay conditions and other regularity 
issues is the language of weighted function spaces; specifically weighted classical 
and Sobolev spaces. 

We will use the definitions of Jl2j , other accounts of the theory may be found in 
0, @ and g§. 

Definition 3. Weighted Sobolev spaces can be defined via the weighted Lebesgue 
spaces Lfi, 1 < p < oo which are spaces of locally measurable functions for which 
the norms 

\\f\\,s=\ ^\f\ P ^ff, P<™ 
[ ess sup^ (a °\f\) , p = oo, 

are finite. We can replace JR n with subsets 51 of JR n in these definitions. The 
weight a is usually taken to be a — yl + r 2 or a = r on subsets excluding {0}. 
The weighted Sobolev spaces are now defined as consisting of functions with weak 
derivatives up to order k for which the following norm is finite 

k 

ll/lk P ,<5 = E W&fWvt-i- 

3=0 

For p < oo we denote the weighted Sobolev space byWg' p . For p — oo we denote 
the classical weighted function space by Cg . 

We will also require the Sobolev inequality, given here in the form presented in 
@- 

Sobolev Inequality 

If F £ Wg' p then 

(i) ||/||^,tf<C||/|k„*. if n~kp>0 and p<q< ^ 
(H) \\f\\oo,S <C\\f\\k, p ,s, ifn~kp<0, and 
\f(x)\ = o(r 5 ) as r — > oo. 
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Another tool we will require is the "multiplication lemma" . In what follows we 
will be mainly using Sobolev spaces with p = 2, and it is for this case that we give 
the multiplication lemma (adapted from Jl3|). 
Multiplication Lemma 

Pointwise multiplication on TR n is a continuous bilinear mapping 
if fei, &2 > k, k < ki + ki — n/2, and S > Si + 62- 

In discussing the decay conditions on ip it will be useful to have some new 
notation. Firstly, we introduce two new 2-spinors oa and t A by, 

ua = VRe^ lx OA 
v A = y/Re^ lx L A: 

(2) with l c o c = 1. 

Here, R is a positive real function and x a rea l function (in fact, as it turns out, 
the same function introduced in the static case). Our non-degeneracy condition 
now reads 

v A u A = Re lx ^ 0. 

The dyad oa , M give a spinor dyad basis which is "co-moving" with the "Dirac 
flow" given by j a . In general we have 

j a = Ro» A (oV + iV), 

with j a j a — 2R 2 . These ideas can be developed further and lead to a 
Newman-Penrose type formalism for the Maxwell-Dirac system (for the 
Newmann- Penrose formalism in General Relativity see Jig]). We will not fully 
pursue this here, although we go somewhat down this path in our proof of lemma 
2, below. 

For our static systems we have 

R = X 2 + Y 2 

V = RV = i?(sin r cos rj , sin r sin 77 , cos r) 
X = \/i?cos — 

Y = Vfisin^ 

( 0A ) = ( sin i e ^; ) 

\ — cos I e^ v J 

(3) w=( sin t e :%) 

\ — cos 2 e 2 1 J 

Notice that oa and l a are both O(l) as r — > 00, so ip decays as VR. 
In discussing decay conditions on the A a we need, of course to be aware that A a 
is defined only up to a gauge transformation. This problem is usually resolved (to 
some extent) by imposing gauge conditions, such as the Lorenz gauge. We restrict 
our attention to stationary gauges and demand that A be 0(l/r) as r — > 00 in 
some gauge. If we demand that A be 0(1/ r) then the question is, can we solve 

V 2 + V.A = 0, 

so that A' = A + is 0(l/r)? That is, can we find a gauge transformation 
which takes us to the Lorenz gauge with the new A still satisfying the appropriate 
0(l/r) decay? The answer is yes provided we choose the original A in the correct 
function space. 
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We need A a to be at least twice (weakly) diffcrentiable to make sense of the 
Maxwell equations. So to get the appropriate differentiability and decay we will 
take (see definition, below) 

A a G W 2 _f +e (E p ), 

for some p > 3/2 and all e > 0. Here, E p = TR 3 \B p , with p < oo large enough so 
that B p encloses all external sources. This means we can now solve the gauge 
equation for <f> G W 3 *{E p ), with A' = A + V(j> G W*f +e (E p ), in fact, for 
< e < 1 the Laplacian gives an isomorphism between the function spaces W 3,p 
and W±f +e , see @. 

The Maxwell equations imply that j a G W^ +e (Ep). In the vector basis 
co-moving with j a (induced from the co-moving dyad) the charge density is V2R, 
we would expect therefore that R is 0(1/ r 3 ) as r — > oo. We also require at least 
three derivatives for R to define the Maxwell equations when the A a are written 
in terms of the components of ip. This suggests we should take R G W_f +e (E p ). 
The oa and i A must also have at least three (weak) derivatives and we also need 
to ensure that j a G Wlf +e (E p ). We will require o A , la G W?> p {E p ), for any e > 0. 
This leaves the differentiability and decay of x to be determined. Again we will 
require at least three derivatives of \- The decay rate, however must be 
determined from the equations. 

We can now give our definition of an isolated Maxwell-Dirac system. For 
concreteness and ease of manipulation we will restrict our attention to Sobolcv 
spaces Wg' 2 (E p ). 

Definition 4. A stationary Maxwell-Dirac system will be said to be isolated if, in 
some stationary gauge, we have 

* = e^R ( Zl* ) , 



with E constant and 

R G Wl^ +e (E p ); oa, la G W 3 ' 2 (E P ) and A a G Wlf +e (E p ), for some p>0 
and any e > 0. 

Remark 

This definition ensures, after use of the Sobolev inequality and the multiplication 
lemma, that ij; = o(r _ 2+ e ) and A a — o(r~ 1+e ). 

We are now in a position to prove a lemma which will be used in the proof of 
theorem 1 in the next section. But we first need some new notation. 
We introduce the complex null tetrad vectors 

r = v a AA <> A d A , n a = or a AA L A Z A 
m a = a a AA o A Z A , rfi a = a a AA i A o A . 



This null tetrad can now be used to define the following (Newman-Penrose) 
intrinsic derivatives 
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With this notation and the expression for A AA of section II we find that the (real) 
potential A a has to have the following components with respect to the null tetrad 



(4) 



A, 
An 

4 

A _ 
i in 



] A a — 

n a A a = i 
2e 

m a A a = i- 
le 



m a A a = — 
2e 



\p2m cos x — A% + i(n — /i + 7 — 7) 
\p2m cos x + Dx + i{p — p + e — e) 



SR 

—g + a-0- 

SR - 
~R- a + p - 



Here a, (3, r, //, p, 7 and e are the NP spin coefficients (Ricci rotation coefficients 
for the non-holonomic NP tetrad), see [fl6|| . Their exact form is not important 
here, what we do need to know is that they are all of the form odo, idt, odi or bdo, 
where d is any one of the NP intrinsic derivatives, D, A, 5 or 6. Notice that here 
we are using a gauge which has the factor e lEt of definition 4 removed - for the 
stationary case this means all variables are independent of t and ^4° — * E/e as 
r — > 00. 

Lemma 2. For an isolated, stationary Maxwell- Dirac system the following must 
hold 



5R/R - |m 
5R/R- frh 

for any e > and some constant E. 

Proof. A straightforward application of the multiplication lemma. 

With o A ,iA& W^ 2 (E p ) and A - E/e, A E W 2 f +e {E p ), for any e > 0, we have 

A m - m — , A m - fh — e W s ' (E p ), 5i > -1 + 3e 
e e 1 

and 

a,/3,7,r,M,7r,e G Wf 2 ' 2 (S^), <5 2 > -1 + 4e. 
The result then follows from (^). 

IV. Vanishing Total Charge 

We will be working with the stationary, static Maxwell-Dirac equations. From 
section II they are 

(Vx)xV = 

V.V = -2m(X 2 + r 2 )sinx 
o m (V X ).V 

= - C0SX+ 2e(Xi + Y*) 



□ 



A = 



1 



[(X 2 - Y 2 )Vn - V X V] 



2e(X 2 + Y 2 ) 
V 2 A° = 47rei° 
(5) Vx(VxA) = 0. 

We can now state and prove our theorem of vanishing total charge. 



Theorem 1. An isolated, stationary, static Maxwell-Dirac system is electrically 
neutral. 
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Proof. A more restricted version of this theorem was proved in fl7|| by one of us 
(H.B.). 

The stationary gauge of definition 4 is the one for which A — > as r — > 0, the 
stationary gauge used in equations (0) has ip independent of t. A gauge 
transformation of the type ip — > e~ ij) will bring the gauge of definition 4 into 
that of equations (||). The A of these equations then differs by a constant, —E/e, 
from the A of definition 4. In fact, we will be able to determine E in the proof of 
the theorem, see also corollary 1. As we noted earlier we have 

OA = 
l a = 



in the static case. The system is isolated so using the multiplication lemma we 
have that 

sin r , cos r , sin -q , cos r\ &W^' 2 (E p ). 
Now we have too = Gq A oa^a = 0, so using lemma 2 we have 

f and f e Wlte(E P ). 



Which in our static case gives, 



, dR . dR\ , dR - a n , 
cost ( cos?? — +sinr? — 1 - smr — G W_' 1+e (E p ) 



(6) and sin rj ^ - cos rj ^ € W!f +e 
The second of equations (Q) can be written as 

(7) -2msinx = = — ^— + V.V, 

with V = (sinr cosry , sin 77 sin r , cost). Using the multiplication lemma we have 

V e Wt 2 (E p ) and V.V G W!f +2e (^), 
for any e > 0. We also have 

V.VR . dR . . dR dR 

(8) — - — = sm t cos 77 — h sin t sm 77 — — h cos t — — . 

R ox ay az 

Next we utilise the invariance of our equations under Lorentz transformations. In 
fact, we know that if (uA(x a ),VA(x a )) is a solution to the Maxwell-Dirac 
equations then (uA{% a ), VA(x a )) is a solution to the original system in the x a 
coordinates; here x a are the Lorentz transformed Cartesian coordinates. This is, 
of course, true for any linear Lorentz invariant theory. Consider the rotation 

x = x cos lu + y sin uj 
y = ~x smuj + y cos a; 

z = z, with u) = . 

2 

This gives 

sinr/(£ a ) dR(x a ) cosry(x Q ) dR(x a ) _ cosr](x a ) dR(x a ) smri(x a ) dR{x a ) 
R(x a ) dx R{x a ) dy _ R(x a ) dx + R(x a ) dy 
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Diffcomorphisms IR 3 — ► IR 3 induce an isomorphism of the Sobolev spaces Wg' p , 
see . The rotation above preserves E p and will give an isomorphism of the 
Sobolev spaces W s ' p (E p ). Using our last equation and (||) gives 

cos 77 dR sin 77 dR 2 2 , . 

This equation and (Q) with the multiplication lemma then gives (multiply each 
equation in turn by sin 77 and cos 77 etc.), 

ldR ldR 2 , 2 
Rte'R-d^ £W -^ E ">- 

The rotation 

x = —z, y = y z = x 

gives 

smri(x a ) dR(x a ) cosr](x a ) dR(x a ) sin77(i Q ) dR(x a ) cosr)(x a ) dR(x a ) 

R(x a ) dx R(x a ) dy ~ R(x a ) d% R(x a ) dy 

Again using the multiplication lemma with (^), we have 

sin 77 dR 2, 2 
~~R~ ~dz — p)' 
In the same fashion, using the rotation x — x, y = z, z = —y, we have 

cost; OR 2i2 

Another use of the multiplication lemma with the last two equations gives 

\dR/dz£W% e {E p ). 

Altogether we have 

ldR ldR ldR 2)2 
Rdx' Rdy' Rdz^ W -i+^r> 
A final use of the multiplication lemma with (^|) and we have 

V.VR ,2 , * 

— ^- G w!f +e (s p ). 

We can now conclude from ([j]) that sin \ G W_'i +e (^ p ), for any e > 0. By the 
Sobolev inequality sin y = o(r~ 1+c ) as r — > 00. 

Now sine is an invertible C°° function on the range of x (on E p , with p large) 
with sin% = for y = 0. So we can now write 

X = tit: + (i, with n = 0, ±1, ±2, . . . and 

(9) Me^W 

Next we use the first of equations (@) to rewrite ^4° entirely in terms of x and 
|Vx|- This equation implies we may write V = 7VX, for some function 7. We 
also have \ V\ = X 2 + Y 2 = \y\ |Vx|, so that 

A = -cosx +— |V X |, with e=^. 
e 2e |7| 

Which, using ([)]), may be written as 

*r> m , s 2to . 2 u £ 

^4 cos(777r) = cos(n7r) sm (— ) + — Vu . 

e e 2 2e 



Hence, 



A - -cos(mr) G V^ 2 2 +£ (£ p ). 
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A° - — cos(mr) G W^ +e (E p ) 



Note that |Vx| is bounded on E p . So 7 cannot change sign on E p as 
R = X 2 + Y 2 ^ 0, e is fixed. 

From the (first) Sobolev inequality applied to ^4° — — cosn7r, with p = q = 2, we 
have that 

A°~ ™cos(ror) G W°_'* +e (E p ). 

Consequently, we also have 

A -— cos(nTr) €W°i(E p ). 
e 4 

Now, in the static case, j° — y/2R, so V 2 (A° — cosn7r ) G W_' 3+e . Hence, 

cos(mr) G W a5 '" 

e 

and from the Sobolev inequality we find that |<9iA°| < Cr~ 2+e (with dt = d/dx 1 ). 
So we have, 

| 9 ^Q|6 r -6(-l-|)-3 < c 4| 9 . A 0|2 r -2(-3+ e) -3 ; 

for any < e < 1/12. Thus, d.A G W"°' 9 6 /4 (£,,). Which finally gives us, 

A a - -cos(mr) G W^'fCE-), 
e _ 4 

and the Sobolev inequality now gives 

A a — cosn7r = o(r~ 3 ). 
From which it is clear that the total electric charge of the system 

lim — / (VA°) .dS, with S p the sphere of radius p, 
p^oo 4tt J Sp y 

must vanish. □ 

Corollary 1. In the gauge in which A — > as r — > 00 the Dirac bi-spinor of an 
isolated, stationary, static Maxwell- Dirac system takes the form 

i{j = e ±im V, with G W^l, (E p ) (as above). 
Proof. This result is a simple consequence of the proof of the theorem. We had 

A -^cos(nTr) G W^ +e (E p ) . 
The constant term, — ^ cos n7r is removed by the gauge transformation 

tb — > e - icos ( n7r ) mt ^,_ 

□ 

In j|] we presented a unique spherically symmetric solution which provides a good 
example of the properties just described. For large r we found 

1 1 , 1 . 

* = 77 ira 3 3 + ° "6") 

mr loom' 3 ?" 3 r D 

TO 1 3 .,1. 

^ = + 9 - 777; w + 

e emr Wlemrr^ r° 

V. The Axially Symmetric Case 

By axially symmetric we mean that the system is invariant under rotations about 
a fixed axis. This requires that gauge invariant quantities be invariant under 
translations in the azimuthal coordinate d>. 
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Definition 5. A Maxwell-Dirac system will be called axially symmetric if 

= 0, 



where 



L = 
N 



d 



„<* k u A u A — 

AA dx a 
AA dx a 



and 



For our static systems we require only that L be invariant under translations in 
In fact, writing L in cylindrical polar coordinates, 

1 'v p p + V+$ + V*k) 



L = 



L dt 



= L° 



d_ 

dt 



V2 

V5\ 



V" 



d 



y<t> 



d 



d 



we find that 



dp 

V = 2XYcos(n- 



= 

must all be independent of 
axially symmetric if 



X 2 + Y 2 

2Xrsin(ry-0) and V z = X 2 - Y 2 
4>. This means our static Maxwell-Dirac system is 



dX _ dY _ d{r, - <j>) 
d(j> dcf) dcj) 



0. 



This information lets us characterise stationary, axially symmetric, static 
Maxwell-Dirac systems as follows. 

Lemma 3. A non-trivial static, axially symmetric Maxwell-Dirac system is 
stationary if and only if rj = <fi in the gauge given in (|l|) . 

Proof. As V is independent of 4> then so is x, as sin y = (V.V r )/(X 2 
reality condition, 

dV 

(V X )XV = 



Y 2 ). The 



dt 



gives, 



~df 



dVP 


- v'f'— = 


dt 


dz 




+ = 


dp 


dz 


dV z 


+ v*^ = 


dt 


dp 





0. 



If the system is stationary lemma 1 says V is independent of t. So either = 
or x is constant. Constant x leads to the trivial solution with X = Y = 0. So we 
must take = 2XYs'm(r] — 4>) = 0. We have rj = <f> mod rm. 
On the other hand, if r] = cj> we have = and consequently both V and V z 
are independent of t. Combining these with dt(X 2 + Y 2 ) = gives the result. □ 

From now on we study the axially symmetric, stationary, static Maxwell-Dirac 
equations. 

It will prove convenient in what follows to use spherical polar coordinates (r, 9, <p) 
and to make the following change of variables 

X = VScos(^) and Y = Vfisin(^). 
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All our dependent variables depend only on (r, 0) , the equations are 
V = R [cos(r - 6) f + sin(r -6)0 

(v x )xv = o 

V.V = — 2mi?sinx 



(10) 



A = — cos x 

e 

2e 1 r sin 



|Vxl 



2e 

cos r 1 
rR 



^-(rRsm(T-e))~^(Rcos(r-9)) 



0- 



together with the Maxwell equations. Here r, and </> are the unit coordinate 
vectors. 

We note that A a automatically satisfies the Lorenz gauge condition. 

In the spherically symmetric case || we have r = 9 with R and x functions of r 

only. In which case 

1 cot 9 f 
2e r 

the magnetic monopole. 

Other tractable cases are those for which r is constant. It is straightforward to 
show there are really only two cases, see p7[ , 

• the cylindrically symmetric case, t = 7r/2, sec J10[. 

• the case r = 0, variables depend on z only, see p8fl. 
VI. Magnetic Monopoles 

The spherically symmetric solution has an external (i.e. not sourced by the Dirac 
field directly) electrically charged magnetic monopole. In this section we will 
prove a theorem which shows that this is, to some extent, the generic situation. 
We will show that the axially symmetric, stationary, static Maxwell-Dirac system 
can have an external Coulomb point charge only if it is magnetically charged. 
First we define what we mean by an external Coulomb field. 

Definition 6. We will say that a Maxwell-Dirac system has an external Coulomb 
field if we can choose spherical polar coordinates and a ball B p centred at r = 
such that 



r 



h, in Bp, p > 0, 
with, h, a bounded function on B p and q constant. 



Remarks 

1. q/r is of course harmonic on B p \{0}, so it is not directly sourced by the 
Dirac field via the Maxwell equation W 2 A° = Airej . In this sense the 
Coulomb field is "external" to the Dirac field. 

2. The condition that h is bounded on B p is quite weak. In practice it will 
follow from elliptic regularity of the Poisson equation V 2 /i = 4irej° 

(j° = V2(X 2 + Y 2 ) = V2R must be at least L l (B p Y for the total Dirac 
charge to be well-defined on B p ), see for example [[l9| and ^0|. We also 
require R to be differentiable at least three times (if only in the weak sense) 
to satisfy the Maxwell equation for A - this puts R in W s ' p for some p > 1 . 
Consequently, h will be in W^f 2 ; which ensures that h can be included in 
one of the classical weighted function spaces. 
Now for our theorem. 



Theorem 2. Suppose an axially symmetric, stationary, static Maxwell-Dirac 
system has an external Coulomb field. Let A( p<pi \ = B P \B P1 , with p > p\ > and 
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Bp as in the definition above. If rd r T — dgR/ R is bounded on Ai PtP1 \, then the 
Coulomb point charge necessarily carries a magnetic charge of Dirac value, ^j- ; 
i.e. all electric point charges also carry a magnetic monopole. 

Remark 

The condition rd r T — dgR/R is bounded on A^ ppi j is not particularly strong. It is 
true, for example, if X, Y € Cq(A( PiPi )) - we assume R = X 2 + Y 2 ^ 0; or, when 
A G Wl 2 (B p \{0}). 

Proof. In jl7j it was shown that if a Maxwell-Dirac system has a central Coulomb 
charge then the magnetic field is necessarily unbounded at the centre. 
In our proof of theorem 1 we noted that V = 7VX, so from ( [To| ) we have 

&rX 

cos(r — 9) = e — 1 and 

(11) sin(r — 8) = e , , where s = ^-r. 

We will first show that sin(r — 9) — > as r — > 0. 
Now, A = £ + h, so from ( |l0|) we have 

(12) \V X \ = g +9, 

r 

where g = h — — cos\ is bounded on B p . Next, we write 

X = q\nr +(, on B p \{0}. 



Equation (112) is 



2 , 2a » a 1 ra s\2 _ 2( ? , 2 



(13) (ft-cr + 7^c + ^(*cr = fs+s 

Note that r|V%| — » <? as r — * so that rd r x = q + rd r ( and dgx = ^6»C are 
bounded on Ai PiP1 \, for p small enough. Let 77 be the smallest non-negative 
number such that r^drC, — > as r — > 0. We can write (|l^) as 

{r^drC) 2 + Zqr^drC + {r^ dgC) 2 = 2qr r 'g + r" +1 g 2 . 

It is clear that r _ 5~ l <9 r .£ — + and r 3 ~i~dg( — > 0. From which we conclude that 
77 < 1 and so dgx = dg( — > as r — > 0. Consequently, 

sin(r — 6) = , — > 0, as r — > 0. 
V ; r|Vx| 

Now rd r x — * q as r — > 0, so r<9 r X cannot change sign as 9 varies, for p small. 
Writing S\ = q/\q\ we have, from ([Tlj), 

„,v cos r(r, it) — > — eei 1 . 

14 V A as r^0. 

v cos r(r, 0) — > eei J 

We now work on A( p , Pl ) with p small enough that cos(r — 9) 7^ 0. The magnetic 
charge of the magnetic field B = V X A in B r , p > r > pi , is 

(15) 6 =7-/ B -dS=\l d e (rsm9A)d9, 



47r is r 2 J 9=0 



where A = A<\> is given in (10). After some manipulation the third equation of 
( |To| ) can be written as 

dgr — 2 H — = tan(r — (?) ^r<9 r r — J — 2mrsinx — 3cos(r — 0). 
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So that if rd r r — is bounded on v4( PlPl ) then so is der — 2 + From (|Io| ) 

we have 



r sin 6 A = — jcos r — P sin } , 
2e 1 ; 

where -P = I ro r r — I cos(t — 0) — I c%t — 2 H 



sin(T - 6>). 



Clearly, under the conditions of the theorem, P is bounded on Ar p<pi ) . From ( |15[) 
we now have 



cosr(r, 7r) — cos r(r, 0) 



cosT(r, 7r) — cos r(r, 0) 



[P sin^] J : 



9=0 



1 

Ye 

Finally, from equations (|l4|) and (|15|) we obtain the magnetic charge in the limit 
r -> 

-££l ±1 



&0 



2e 2e 



□ 



Corollary 2. Suppose we have an axially symmetric, isolated, stationary, static 
Maxwell- Dirac system with the only external sources being N isolated electrically 
charged magnetic monopoles. Let the conditions of theorem 2 apply in the N balls 
B Pi , containing the charges. Let the conditions of theorem 1 apply on Ep . Then, 
if N is even there are N/2 positive charges and N/2 negative charges (with 
corresponding monopoles), with the total magnetic charge of the system being zero. 
Lf N is odd there are (N — l)/2 charges with one sign and [N + l)/2 charges with 
the opposite sign, with the total magnetic charge of the system being ±l/(2e). 

Proof. We have N charged monopoles each in a ball B Pi , z = 1,2,3 ... N, there 
are no other external sources and all the B Pi are properly contained in B p . 
We have A is 0(l/r) so that P, as defined in the above proof, is bounded on Ep). 
Using the results in the proof of theorem 2 and the divergence theorem we find 
the total magnetic charge of the system. 



If 1 

&tntal = lim ~T- / ( V x A ) - dS = 7T lim 

total r _ >00 4^ y i 2e r-^oo 

N 
i—1 



cos t (r, 7r) — cos r(r, 0) 



This gives, 



N 

1 > £^ £l - _1 ' 



i=l 

from which the results of the corollary follow. □ 
VII. Conclusions 

In classical physics one expects stationary or static systems to be the end point of 
some time evolution. This clearly cannot be the case for a single isolated electron 
modelled by the Maxwell-Dirac system. To construct such a model we will have 
to abandon one, or both, of the static and stationary assumptions. The electron 
must be a dynamic object. The same also applies to other stable objects such as 
the Hydrogen atom. A Maxwell-Dirac model of such objects must be non-static or 
non-stationary, or both. 
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Theorem 1 is remarkable in many ways. No matter what arrangement of external 
electric and magnetic fields inside the ball B p , no matter what we do to the 
topology in B p the total electric charge of the system must be zero. The total 
charge vanishes purely as a result of the asymptotic decay and regularity 
conditions. 
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